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The zero momentum sectors in effective theories of QCD coupled to pseudoreal (two colors) 
and real (adjoint) quarks have alternative descriptions in terms of chiral orthogonal and symplectic 
ensembles of random matrices. Using this correspondence, we compute correlation functions of Dirac 
operator eigenvalues within a sector with an arbitrary topological charge in a presence of finite quark 
masses of the order of the smallest Dirac eigenvalue. These novel correlation functions, expressed 
in terms of Pfafhans, interpolate between known results for the chiral and quenched limits as quark 
masses vary. 
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I. INTRODUCTION 

Random matrix theory of disordered Hamiltonians 
jl],D relies upon an ansatz that in the ergodic regime 
where the system size L is much larger than the elastic 
mean free path but much smaller than the localization 
length, details of the Hamiltonian are lost except for its 
time-reversal and internal symmetries. This ansatz was 
materialized by Efetov H who derived Wigner-Dyson 
statistics out of Anderson tight-binding model by retain- 
ing only the zero-momentum mode of his spectral non- 
linear a model (NLcrM). 

By taking the flavor symmetry among quarks into 
consideration as an additional internal symmetry, Ver- 
baarschot and collaborators @-^| have reinterpreted this 
wisdom in the context of quantum chromodynamics 
(QCD). In this context, the spectral NLcrM is a super- 
symmetric extension [j^U of the conventional NLcrM over 
the coset manifold associated with the chiral symmetry 
breaking 

Z (2) (6>; Af) = / dU exp(Re tr e l6/N f MU r ), (la) 
Jsu{N f ) 

Z {1) {6;M)= [ dU exp(Re tre l9/N f MU3U T /2), (lb) 

J SU{2N f )/Sp{2N f ) 

ZW{$;M)= [ dU exp(Re tr e 19 ^^^ MUU T ), (lc) 

JsU(N f )/SO(N f ) 

after retaining only the zero mode. Here the superscripts 
(2, 1, 4) of Dyson indices (3 refer to the anti-unitary sym- 
metry of Euclidean Dirac operators (which are considered 
as stochastic Hamiltonians) in three classes of QCD [ pd| : 

(3 — 2 : N c > 3, Nf fundamental Dirac fermions, 

(3 = 1 : N c = 2, Nf fundamental Dirac fermions, (2) 

(3 = 4 : N c > 2, Nf adjoint Majorana fermions. 



The rescaled quark mass matrices M are defined as 

M= <iiag(jMi,...,n Nt ) (/3 = 2,4), 

M = diag (/ii, . . . , fi Nf ) ® J (/3 = 1), (3) 

Hi = Y>L A mi, J = l Nf (g) J, ^ = ^ ) ' 
with a limit 

L — > oo, mi — > 0, fii : fixed (4) 

being assumed, £ stands for the quark condensate in 
the chiral limit, and 9 stands for the vacuum angle. 
By the same token as the spectral NLcrM of the tight- 
binding model was derived from a conventional ran- 
dom matrix ensembles |l2| , these NLctMs have an alter- 
native derivation from chiral random matrix ensembles 
(xRMEs) OH: 

zC«({m})= r 'f e -'"W(0;{m}) 

JO 27T 

Nf 

= J dWe~^ v ^ n det i^ZX (5) 

i— 1 ^ ^ 

under a limit 

N —>■ oo, rrii —>■ 0, [ii — Trp(0)rrii : fixed. (6) 

Here the integrals are over complex, real, and quaternion 
real (N + v) x N matrices W for f3 — 2, 1, 4, respectively, 
and p(0) stands for the large- N spectral density of the 

random matrix T> = ( ^Q: 

p(X) = lim (txS(X-W)), (7) 

N—>oc 

at the origin. It is understood for f3 = 4 that twofold 
degenerated eigenvalues in the determinant arc only 
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counted once, and the topological charge v is substi- 
tuted by vN c . These xRMEs are motivated by the mi- 
croscopic theories (Euclidean QCD) on a lattice, with 
a crude simplification of replacing matrix elements of 
the anti-Hcrmitian Dirac operator^) = [d^ + iA^)^^ by 
random numbers T> generated according to the weight 
e -/5tr v(w 1 w)_ u nc Jer this correspondence, the micro- 
scopic limit (^|) is equivalent to Leutwyler-Smilga limit 
(f|), since the size N of the matrix W is interpreted as 
the number of cites L 4 of the lattice on which QCD is dis- 
cretized, and the Dirac spectral density at zero virtuality 
p(0) is related to the quark condensate by Banks-Casher 
relation £ = np(0)/L 4 JIJ]. Using the xRME represen- 
tation, various correlation functions of microscopically 
rescaled Dirac eigenvalues, 



N — > oo, A -> 0, C = 7rp(0)A : fixed, 



(8) 



have been computed for all three values of (3 in the mass- 
less case fi = [p~5| p3| . On the other hand, in a pres- 
ence of finite /x's, Dirac eigenvalues comparable to fi's 
are expected to obey statistics that interpolate the chiral 
(/i — ► 0) and quenched limits (/i — > oo or Nf — 0). These 
novel spectral correlation functions have been analyti- 
cally computed, until recently, solely for the chiral uni- 
tary {(3 — 2) ensemble |24| |2fJ. Therefore we aim to treat 
the remaining cases, chiral orthogonal (/3 — 1) and sym- 
plectic ((3 = 4) ensembles with finite mass parameters, 
and compute Dirac eigenvalue correlation functions for 
these ensembles. We anticipate that advances in numer- 
ical simulations of lattice QCD with dynamical quarks 
p5|,pT| will confirm our analytical results, Eqs.(|48|) and 
( p3| ) below, in a foreseeable future. 

This Article is organized as follows. In Sect. 2 we 
compute the correlation functions for the chiral orthog- 
onal ensemble, by utilizing the quaternion determinant 
method developed in Ref. |Q. In Sect. 3 we exhibit an 
explicit form of the correlation functions for the chiral 
symplectic ensemble, which was obtained by the Authors 
p3[ as a corollary to the computation of the partition 
function. In Appendix A we collect definitions related to 
a quaternion determinant. In Appendices B and C we 
present alternative expressions of the quaternion kernels 
for the orthogonal ensemble, and the symplectic ensem- 
bles with quadruply degenerated masses, respectively. It 
enables us to identify our results with those in a paper 
by Akemann and Kanzieper |$4j, which appeared after 
the Letter by the present Authors |33l] and computed a 
1-level correlator with a single massffor [3=1) and p- 
level correlators with quadruply degenerate masses (for 
= 4). 



II. ORTHOGONAL ENSEMBLE 

We start by expressing the partition function (|5|) of the 
XRME in terms of eigenvalues Xi = Af of the Wishart 
matrix W^W (up to an overall constant): 



Z { f ) ({m})= (n™r)So({m}) 



S ({m}) 



oo />oo 



N 



JO J = 1 
N 



Y\d X] P({x}-{m}), 



(9) 



P({x};{m}) \[L ' V ' \[ir, • nr 

j=l \ i=i 

N 

x Y[ \Xj-X k f. 

j>k 



(10) 



The indices (3 and v are suppressed for simplicity. Since 
the partition function (^) is even under v i— ► — u, we have 
set v non-negative integer, without loss of generality. The 
p- level correlation function of the Wishart matrix W^W 
is defined as 

a(xi,...,x p ;{m}) = — — - , r _^ , (11) 



5 ({m}) 



E p (xi,...,x p ;{m}) = 

/•OO />CO N 



-i pOO pOO 

— ^ / ... / [] d Xj P({x};{m}). (12) 

-py-Jo Jo ]=p+1 



(N- 



Then the p-level correlation function of the block off- 
diagonal Hcrmitian matrix iT>, 



p(Ai 



,\ P ;{m}) = (i[tv5(X k -iV)), 



(13) 



k=l 



is expressed in terms of a multiplied by the Jacobian of 
the transformation A i— > x — A 2 : 

p(Xi, . . . , A p ; {m}) = 2? JJ |A,| a(\j, . . . , A 2 ; {m}). (14) 

3=1 

As the universality of correlation functions of the uni- 
tary ensemble in the microscopic limit (||) |35| [$7j,|o| is 
known to inherit to those of orthogonal and symplectic 
ensembles 38 13^], it suffices to concentrate on Laguerre 
ensembles, V(x) = x. This leads to Wigner's semi-circle 
law 



p(A) = —\J2N — A 2 



(15) 



Now we concentrate on (3 = 1 , and define new variables 



. a, 



-m 2 (<0), j = l,. 
Xj-a (> 0), j = a + 1, . . . , a + N. 

Then the multiple integral (|lj) is expressed as 
"p(^l, • • ■ , z a+p ) = —— - — — ■ - 

roo r oo a+N a+N , a+N 



(16) 



x (17) 



(N 



| p OO p oo _ i 

r^r - 1 n ^nv^n 

yj jo jo J=a+p+l j= i j>k 



Zj - z k , 
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where 



The quaternion elements f(z, z') are represented as 



w (z) = |«| tf-1 e- 



(18) 



Eq.(|l7|) resembles an (a + p)-\eve\ correlation function of 
the conventional (massless) Laguerre ensemble with a + 
N levels. However, conventionally the levels z\,..., z a+p 
are all positive, while in the present case some of them 
(zi, . . . , z a ) are negative. We carefully incorporate this 
fact into the following evaluation. 

Let us denote the integrand in Eq.([l7|) as 



7 



P&l, 



7 



' z ^ = n v w ( z ^ n i z i ~ zk 

j=l j>k 



(19) 



with 7 = 0; + ^. It can be readily seen that 

7 7 

p{zi,...,z 1 ) = JJ ^Jw(zj) Y[{zj - z k ) x 

3=1 j>k 

Pf[sgn(z fc - Zj)] jtk=li _„ 

[sgn(z fe - Zj)}j,k=i,... n [ft]i=l,..., 7 
[-.9fc] fe =i,..., 7 



Pf 



(20) 

(7 : even) 
(7 : odd), 



with gj = gk = 1. The Pfafnans in the above can be rep- 
resented as quaternion determinants p0|- p2| , p2|p^j32] , |43[ | . 
In doing so, we need to introduce monic skew-orthogonal 
polynomials R n (z) = z n + ■ ■ •, which satisfy the skew- 
orthogonality relation: 

{Rzni R-2-m+l) R = — (R2rn+1 , i?2n) R = f n 5 nm , Others = 0, 

(21) 

where 



(f,g) R = dz^Mz)g(z) dz'^Mz^f(z')-(f^g), 
Jo Jo 

(22) 



and (integrated) 'wave functions' 



^n(z) = y / w(z)R n (z), 



(23a) 
(23b) 

(24) 



$„(z)= / dz'sgn(z - z')^w(z')Rn{z'). 
Jo 

Note that for negative z, & n (z) is a constant: 
$„(z < 0) - $„(0) = -s„. 
Now we present the following theorems: 
Theorem 1 

For even 7, we can express p{z\, . . . , z 7 ) as 
,7/2-1 

P(*i, •••)%)= ( J! r j ) Tdet [/( 2 i. z fe)kfe=i,...,7- (25) 



f(z,z') 



S(z,z') I(z,z') 
D(z,z') S(z',z) 



(26) 



The functions S(z, z'), D(z, z') and I{z, z') are given by 

, (27a) 

.(27b) 

(27c) 



7/2-1 

S(z,z') = E 
n=0 

7/2-1 

D(z,z') = ]T 

n=0 

7/2-1 

I(z,z') = J2 

n=0 



$2n(z)*2n^ 


-i{z')-® 2n+1 {z] 


*2„(*') 




r„ 




L *2n(-Z)*2n- 


hl (z')-*2n+l(^ 


)*2n(*') 




r„ 






-l(z')-*2n+l(2) 


$2„(^') 





Theorem 2 

For odd 7, we can express p(zi, . . . , z 7 ) as 
p(zi, ...,z 7 ) = 

/[7/2]-l v 

n r J -)a 7 -iTdefc[/ odd (z i ,z*)k*=i 
The quaternion elements are represented as 
/° dd (z,z') = 



(28) 



S^^z') 7 odd (z,z') 
D odd (z,z') S odd {z',z) 



(29) 



and s„ is defined in Eq.(§|). The functions S odd , L> odd , 
and / odd are given in terms of S, D, and I in Theorem 1 
according to 



S odd {z,z') =S(z,z') 
D odd (z,z')=D(z,z') 
I° dd (z,z')=I(z,z') 



*7-l(Q , 

S7-1 



$ 7 _i(z) - $ 7 -i(z') 



S 7 _l 



Here * stands for a substitution 

s 



R n (z) i-> i?„(z) — Ay_l(^) 

s 7 _i 



(30a) 
(30b) 
(30c) 

(31) 



for n = 0, . . . , 7— 2, associated with a change in the upper 
limit of the sum 



7/2 - 1 ~ [ 7 /2] - 1. 



(32) 



Theorem 3 



Let the quaternion elements qjk of a selfdual nxn matrix 
Q n depend on n real or complex variables z±,---,z n as 



Qjk = f(Zj,Z k ). 



(33) 
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We assume that f(z, z') satisfies the following conditions. 

f(z,z)dfi(z) = c, (34) 
f(z, z")f(z", z')d t i(z") = f(z, z') + Xf(z, z')-f(z, z')X. 



Here d/i(z) is a suitable measure, c is a constant scalar, 
and A is a constant quaternion. Then we have 



/ 



Tdet Q n dfi(z„) = (c - n + l)Tdet Q„_i, (35) 



where Q n -i is the (n — 1) x (n — 1) matrix obtained by 
removing the row and the column which contain z n . It 
is straightforward to show that the quaternion element 
/(z, z') and / odd (z, z') in Theorem 1 and Theorem Uboth 
satisfy the conditions imposed on /(z, z') in Theorem 3 
with dfi(z) = w(z)dz. This means that we can write 



E p (zi, 



) = ■ 



n [( a +N)/2]-l, 
llj=0 



■x (36) 



117=1 n°> fe i z i - 

Tdet[/(z 3 ,z fc )] :) -. fc= i.... !Q+ p (a + N : even) 
s^-iTdetl/ ^^-,^)]^^!,...,^ (a + iV:odd) ' 

Since the final result in the asymptotic limit N — > oo 
should be insensitive to the parity of N, we consider only 
even a + N henceforth. Then the p-level correlation func- 
tion (111]) is finally written as 



a{xi, . . .,x p ;mi, . . . ,m a ) = 



Tdet[f(zj,z k )]j,k=i,..., a 



Tdct[f(Zj, Z k )]j,k=l,...,a 



(37) 



Now we proceed to take the asymptotic limit of the 
correlation function, by making use of explicit forms for 
the skew-orthogonal polynomials and their norms associ- 
ated with the weight (|l8|) obtained by Nagao and Wadati 



R2n+l(z) 



_ (2n)! d T ( V -1) (2 x 

~ 2 2n + 1 dz 2n+l[ h 

_ (2n +!)! („_!) 

_ 2 2n+1 2,1+1 ^ ' 



(38) 



(2n)\ 



^(2n + ,)-4r 1) (2.) 



2 2 



r n = 2 



—Aji — v 



dz 

(2n)\{2n + v)\ , 



expressed in terms of the Laguerre polynomials 



z e 



! dz- 



-(e- z z n+a ). 



(39) 



We need to evaluate the local asymptotics of the quater- 
nion function /(z,z'), whose constituent Laguerre poly- 
nomials tend to 

( } t (n/|.|)° /2 7„(2V^N) (*<0) ' (40) 

as n — > oo, z — > 0, with nz : fixed. Three cases should be 
considered separately: 

(a) z, z' > 0, (b) z < 0, z' > 0, (c) z, z' < 

(the case z > 0, z' < is unnecessary because of the 
sclfduality f(z,z') = f(z',z)). 

(a) z, z' > 

We define microscopic variables £ and (' as 



8JV' 



c' 2 

87V' 



(41) 



according to Eq.(§) with 7rp(0) = 2\/2~/V. If both argu- 
ments are positive, the asymptotic limit is known to be 




S ++ (C,C') = ^S(z,z') J ds fsJu-xW^Qp- J v (ts)J v ^)j + 



■UK') JM) 



(42a) 
(42b) 



(87V) 2 --' -' !<i./,;. v— C C 

7 ++ (C,C') = / d« 2 / / dv(uJ u ^i{tu)J u {tv)-J v {tu)vJ v - 1 {tv))-j J v (u)du - sgn(C - C')- 

■/o Jo Jo JC 

(42c) 



(6) z < 0,z' > 

We define microscopic variables /z and ( as 



'8JV' 



8JY' 



(43) 
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according to Eqs.(||) and (§]). When one of the arguments is negative, the identity (|2J) should be taken into consid- 
eration. We find 



S- + fa, = S + (0 = ±S(z, z') ~ 



£>-+(/*> 



ds I s Ji/_i(ts) — ^ — J u (ts)I v -x(tfx) ) + ' ' 



4/x 



1 



D(z,z') ~ — / ofct 2 ( I v -i{tn) 



c 



J v (u)du + 1. 



(44a) 
(44b) 
(44c) 
(44d) 



(c) 2,2' < 

We define microscopic variables /Lt and Q as 



We can readily derive 

S—(jj,,h') = S-{n') 



£_ 

8N' 



JL S (z zT\ - W) 
8N b[Z ' Z) V : 



87V' 



^ — (Mj M') = J(z> z ') ~ s 8' n (M ~ M')- 



(45) 



(46a) 

(46b) 
(46c) 



In Eqs.(44a), ( |44dD , and (46a), we have introduced symbols with one sign subscript (e.g. 5+) in order to indicate 

that they depend only on the second arguments of those with two sign subscripts (e.g. S |_). 

The scaled correlation function p s is defined as 



> Cpi ■ ■ ■ ? Ma) 



1 / Ci Cp Mi 

"Trvn^v TSttf' * * * ' -7==; -7=) ■ 



(8iV) p V^V' WSNWSN 
By making use of Eqs.(|l4|), (^7|), and Dyson's equality (A6), we finally obtain 

-J_+ 5 

Pf 



'8N 



(47) 



Ps (Ci, . . . , C p ;mi, • • • , m«)= (-l)^- 1 ^ n |ai 



D. 



+ 



fe=i 



Pf 



-S T _ 



= (-i) p(p - i)/2 2 p niai 



k=l 



= (-iy (p - i)/2 2 p i[\(k\ 



k=i 



Pf 




-i+- 

_oT 


- D-+~ 

- s++ 

- D ++. 






Pf [£>__] 













S- 


-1+ 




s+ 


Pf 


-S T 


D— 


sl- 




D-+ 




5+_ 






S++ 




.-4 








D++ 


Pf 










(a : even) 



(48a) 



(a : odd). 



(48b) 



The elements of the matrices S ee i, D ee r, I ee i and the row vectors S e , I c (e, e' = +, — ) in the above are defined as 



5 



(S ++ ) M = S ++ (( k ,Q), (S-+)u = S-+(jh,Q), {S—) i j = S—(jH,fJ'j), (S+)e = S+(Q), (S-)j = S-(fij), etc., (49) 

where the subscripts take their values in i, j — 1, . . . , a and fc, t = 1, . . . , p. In the last two lines we have exploited a 
Pfaffian identity that holds for antisymmetric matrices A, B, and a row vector v: 







V 


Pf 


A 










V 




—v 1 ■ ■ ■ — V 1 


B 



= Pf [A] Pi[B] 







1 






Pf 


A 


1 


Pf 


v 
~v T B 




-1 1 










(rank(A), rank(S) : even) 



(rank(.A), rank(B) : odd). 



(50a) 



(50b) 



In a special case p = a = 1, the expression i 
B). 



reduces to Akemann and Kanzieper's recent result p4| (see Appendix 



-*4 

-SI' D 



+ + u ++ 



In the quenched limit [i\ , . . . , /i Q — > oo when the ratio of two Pfaffians is replaced by a minor Pf 

the correlation function tends to that of Laguerre orthogonal ensemble computed by Nagao and Forrester, Eqs.(2.21 
(2.18), (2.19), and (3.20) in Ref. [p2[ , with v = 2a + 1. By the same token, it satisfies a sequence 



,Ma) 



Ps({(}; Mi 



■ Ma-i) AlQ ^°° p s ({C}; mi, 



, Ma-2) 



(51) 



as each of the masses are decoupled by sending to infinity. To illustrate this decoupling, we exhibit in Fig.l a plot of 
the spectral density p s (C; A*) (u = 0,p = 1, a = 1) that interpolates between known results for the chiral and quenched 
limits. 

III. SYMPLECTIC ENSEMBLE 

Although correlation functions of the massive chiral symplectic ensemble have previously been computed by the 
Authors [^3), we nevertheless present their explicit expressions for the sake of completeness. We concentrate on the 
case with an even Nf(= 2a) number of flavors and pairwise degenerated mass parameters, corresponding to adjoint 
Dirac fermions in the QCD context. The scaled p-level correlation functions, defined in Eq.(|4^), is expressed by 
construction as a ratio of partition functions with 2a and 2a + Ap flavors . 



Pa(d, • • • , tp-, m) = c$i nee* - ci) 4 n (ic*i 3 n^ 2 



fe=i 



2 2\2\ ^ (Ml , Ml i ■ ■ • j Ma , Men iQl, ■ ■ ■ , *Cl > • • • ! *Cpi ■ ■ ■ ? *Cp) 

Mi 



Zv (pi, 111, ... , fX a , Ma) 



(52) 



Here d P i stands for a constant to be fixed below. Using an explicit form of the partition function, Eq.(31) of Ref. 
p3f, and taking confluent limits in Cfc's, we obtain 



Ps (Ci,...,C P ;M)-(-i) p(p+1)/2 2 p ni^l 





J__ 


I-+ 


s. 


-+ 


Pf 


-I T + 


I++ 




-+ 










-+ _ 



fc=l 



pf[/__] 



Pf 





0- 


I-+ 


S-+ 







-Ql 


-Pi 




Q+ 


i++ 


s ++ 




^+ 




D++ 



k=l 



Pf 



/_ Q- 
-Q T 



(a : even) 



(53a) 



(a : odd). 



(53b) 



G 



The elements of the matrices S ee ', D ee i, I cc i and the column vectors Q e , P e (e, e' — +, — ) in the above are defined as 

(I — )ij = I — (Mi) H) = Vi^j / dtt / du(l2u(2tHi)l2u(2tu^j) - I 2v {2tu^i)l 2v {2tnj)), 

Jo Jo 

(/_+)«= 7_ + (/Xi,^) = / dtt I du(l 2v (2tfj,i)J 2u (2tuQ) - l 2v {2tu\i i )J 2v {2tC ! t)) 1 

Jo Jo 

(S-+)ie = S- + (ni,Q) = tM / dtt 2 du(l 2v (2tpn)u J 2u+1 (2tu( e ) - I 2v (2tupn)J 2v+1 (2tQ)), 

Jo Jo 

(/++)« = I++{Ck, <U) = CkCt [ dtt f du(J 2u {2K k )J 2v {2tuQ) - J 2 „(2tu( k )J 2 „(2tQ)), (54) 

Jo Jo 

S++(( k ,Q) = (k / dtt 2 / du(J 2v {2t(, k )uJ 2v+x {2tuQ)- J 2 „(2tu( k )J 2 „ +1 (2t( e )), 
Jo Jo 

(D ++ ) ke = D ++ (( k ,(e)= f dtt 3 [ duu(J 2v+1 (2tCk)J2 V+ i{2tuC i )-J 2v+1 (2tuCk)J2 V +i(2tCe)), 

Jo Jo 

{Q-)j = Mj / dtl 2v {2tnj), (Q+)e = Q f dtJ 2u {2tQ), (P+)t=[ dtt J 2 „ +1 {2t( e ), 
Jo Jo Jo 



where the subscripts take their values in i, j — 1, . ..,a 
and k,£ = I,..., p. The overall constant is determined 
as the above by requiring that in the quenched limit 
fix , . . . , fi a — > oo when the ratio of two Pfaffians is re- 



placed by a minor Pf 



'++ °++ 



the correlation func- 



tion tends to that of Laguerre symplectic ensemble com- 
puted by Nagao and Forrester, Eqs.(2.27), (2.25), and 
(4.7~9) in Ref. |p2[ (whose notations are related to ours 
via an unfolding change of variables 



i++(C,C)-- 
s++(C,C') = 



c 2 c 

-hi—,—), 
n 8A' 8A ; ' 

87V 4V 8iV' 8N h 

i c 2 c 2 > 

J8N^ D ^8N' 8N> 



(55) 



and v = 2a — 1/2). It is easy to confirm that the cor- 
relation functions satisfies the decoupling sequence j5l| ) 
as each of the masses are sent to infinity. To illustrate 
this decoupling, we exhibit in Fig. 2 a plot of the spectral 
density p 8 {Q] jtt, fJ.) {v = 0,p = 1, a = 1) that interpolates 
between known results for the chiral and quenched limits. 

After the Authors announced the above formula in [B3| , 
Akemann and Kanzieper p4| presented another form of 
the asymptotic correlations in a special case of quadruply 
degenerate masses. In Appendix C, we shall reproduce 
their result as a confluent limit of our formula. 
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APPENDIX A: QUATERNION DETERMINANT 



A quaternion is defined as a linear combination of four 
basic units {1, e±, e 2 , e 3 }: 



<7 = <7o + q- e = <7o+ 9i ei + q 2 e 2 + q 3 e 3 . 



(Al) 



Here the coefficients qo,Qi, q 2 and q 3 are real or complex 
numbers. The first part go is called the scalar part of 
q. The quaternion basic units satisfy the multiplication 
laws 



1-1 = 1, 1 ■ ej = ej ■ 1 = ej, j = 1,2,3, 
e i = e 2 — e 3 — e i e 2e3 = — 1- 



(A2) 



The multiplication is associative and in general not com- 
mutative. The dual q of a quaternion q is defined as 



Q = Qo - q • e. 



(A3) 



For a selfdual N x N matrix Q with quaternion elements 
qjk has a dual matrix Q = [qkj]- The quaternion units 
can be represented as 2 x 2 matrices 



e 2 



1 0" 






' 


-1 " 




1 


, ei -> 


1 










—i 






i 


" 


—i 





, e 3 







—i 



(A4) 



We define a quaternion determinant Tdet of a selfdual Q 
(i.e., Q — Q) as 



TdetQ = ^(-l) JV -' l[(q ab q bc 



(A5) 
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where P denotes any permutation of the indices 
(1,2,..., TV) consisting of I exclusive cycles of the form 
(a — > b — > c — » ■ ■ ■ — > d — > a) and (— l)^ - ' is the parity 
of P. The subscript means that the scalar part of the 
product is taken over each cycle. Note that a quaternion 
determinant of a selfdual quaternion matrix is always a 
scalar. The quaternion determinant can as well be rep- 
resented by the 2N x 27V representation C(Q) |io|: 



TdetQ = Pf[JC(Q)], J = l N (g> 



1 

-1 



(A6) 



S(z,z') ~ 2iV 



2 / dttJ v ^(tOJ v -i{tC') (Bl) 
o 

ds J v -2(s) - 1 

o 



Ju(C 
C 



where we have adopted the microscopic variables (|41 
This asymptotic relation was derived by Forrester, Na- 
gao, and Honner |43|| in a study of parametric random 
matrix ensembles. Note that the first integral in the 
above is equal to the Bessel kernel, 



APPENDIX B: QUATERNION KERNEL FOR 
THE ORTHOGONAL ENSEMBLE 

Recently Akemann and Kanzieper 0] derived the 
asymptotic correlation function in a special case p = a = 
1. Though their result is clearly in agreement with that 
special case of ours, there is a difference in the appear- 
ance because they adopted alternative asymptotic formu- 
las. Their formulas are based on an asymptotic relation 
for z, z' > and v ^ 1 (see below): 



CJ,(C)J,-i(C , )-^-i(C)C%(C / ) 
c 2 - C' 2 



(B2) 



From the derivation the equivalence of (42a) and (Bl) 
was well established. However, it is worth directly prov- 
ing it here for an unambiguous identification. 

Using the Bessel function identities, we can readily see 
that 



tf dssJ v - 1 (ts)=CMtC) + (y-l)[ dsJ u (ts), 
Jo Jo 

f c f c ds 

t dsJ v (ts) = -J„- 1 (tC,) + J v - 1 (p) + (i/-l) / —J v -x{ts). 
Jo Jo s 



Substitution of Eq.OBJJ) into (42a) yields 



S{z,z') ~ 2N 



£ / dttJ u (tOMK')+ f dtt(j v . 1 (tO-Ju-i(o))J v -i(K') + ^^ 

C Jo Jo s 



v - 1 f 1 f 1 ds 

—j-l dtt ds J u (ts)J u (t£) - (y - 1) / dtt — J 1/ _ 1 (ts)J v _i(iC') 
C Jo Jo Jo Jo s 



By partial integrations, we find 

C f 1 



and 



,, dttJ v {tQJv{K')= I dttJ^itOJu-iitO-Jv-iiO^P 

s Jo Jo s 



i ['dtt [ C dsJ u (ts)J u (tC)= ['dtt [ C — J v _ x {t8)J v -t{1^)- [ -J,_i(s)^p. 
C Jo Jo Jo Jo s J s C 



We substitute (B5) and (B6) into (EM) and obtain 
S(z,z') ~ 

27V 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



2/ dttJ v -i{tC)Ju-i(tC)-Ju-i(0)f dttJ v -!(tC')- ( J u -i(<0 + (v - 1) I — J, 
o Jo \ Jo s 



c 



Again by a partial integration, we find 



c ds 



Ju-i{0 + {v-l) —J v -x{s) = Ju- 



/,-i(o)+ / C « 

Jo 



ds J w _2(s) 



(B8) 
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and thus arrive at the desired result (Bl) provided that v is not equal to 1. Since J^-i(O) = 6 Vt %, the case v = 1 is 
exceptional. 

Similarly, for z < 0, z' > and v ^ 1, an alternative expression 



S{z',z) ~ 27V 



2 / dttJ^itOh-iitf*)-^- 
o 



in terms of the microscopic variables (|43| ) is available. 



APPENDIX C: QUATERNION KERNEL FOR 
THE SYMPLECTIC ENSEMBLE 

In previous ||| and this works, the Authors evaluated 
Dirac eigenvalue correlation functions for the symplectic 
ensemble with doubly degenerate masses. If masses are 
quadruply degenerate, the evaluation of the correlation 
functions is easier because the conventional integration 
method for the massless Laguerre ensemble works with- 
out any modification. In that case, the multiple integral 
we need to evaluate is 



E p (xi,...,x p ;{m}) 
1 



(N-p)Uo 



N 



N 



n ^il 

j= P +i j=i 



(CI) 



N a 



N 



xnnfo+ m <) 4 ni a *- a * i 4 > 

j=l»=l j>k 

where we set Nf = Act and 

4 4 
{m} = ('TOl, ... ,TOl, . . . , TO a , 

The conventional 'massless' theory 
that the correlation functions are written in terms of 
quaternion determinants: 



S2] tells us 



a(x 1 ,...,x p ;{m}) 

ao(«i, 

_ Tdet[/(zj,z fc )]j,fc=i,,.. iCt+ p 
Tdet[/(z J ,z fc )] :) - :fc= i : ... iQ 



(C2) 



l 

1 



S(z,z') S_+(/i,C), 



A 4 



ds J y _2(s) - 1 



(B9) 



where we have adopted the notation (|16|). 

The quaternion function /(z, z') is represented as 



S(z,z') I(z,z') 
D(z,z') S{z',z) 



(C3) 



We evaluate the asymptotic limit of the quaternion func- 
tion f(z, z') in each of the three cases 

(a) z, z' > 0, (6) z < 0, z' > 0, (c) z, z' < 0, 

as in Sect. 2. in terms of the microscopic variables 
CO,/j' defined by Eqs.@, ©, ©, respectively. 

(a) z,z' > 

In the case z, z' > 0, Nagao and Forrester [p2l derived 
the asymptotic limit (|55|), that is 



—S(z,z') S ++ {(,('), 



1 



(8N) 



T D(z,z')~D ++ ({,C'), 



(C4) 



i(z,z')~-i ++ (<;,£), 

where S ++ , D ++ , and 7 ++ are defined in Eq. (j54|) . 
(6) z < 0,z' > 

Using the asymptotic formula for the Bessel function 
((40|), we can similarly treat negative argument cases to 
obtain 



S(z',z)~C / dtt 2 / du(j 2u {2tuQI 2v+1 (2tfi) - J 2u (2tQuI 2v+1 (2tufi)), 
Jo Jo 



1 



(8N)- 



;D(Z,Z') 



f dtt 3 f duu(l 2v+1 {2tufi)J 2u+1 {2tQ)-I 2v+1 {2t^)J 2v+1 {2tuC)), 
Jo Jo 



(C5) 



where <S |_ and 7 |_ are defined in Eq.(|54 

(c) z, z' < 







D(z,z')~- / dtt 3 / duu(l 2v+1 (2tufi)I 2l/+1 {2tn') - I 2lJ+1 (2tn)I 2 „ +1 (2tun')), 



my jo jo 

I(z,z') 1 — (/i,//), 



(C6) 



where J++ is defined in Eq.(|54|). 

We can use Dyson's equality ( |A6| ) to see that the 
above quaternion determinant expression is identical to 
the limit of quadruple mass degeneracy of the general 
formula ( |53a| ) employing Pfafnans. In this case, yet an- 
other equivalent asymptotic formula was recently pre- 
sented by Akemann and Kanzieper ]34j] . Now we shall 
directly demonstrate the equivalence. We should firstly 
note that, under the change of the quaternion elements 

S{z, z') - S(z, z') ee S(z, z') - ^ 



I(z,z')^I(z,z') 



S(z,z")dz", 



d 



D(z,z')-+ D(z,z') = —S(z,z'), 
oz 



where 



W(z) = \z\ v+1 > 2 e~ 2z . 



(C7) 



(C8) 



the quaternion determinant is unchanged. This transfor- 
mation was introduced in Ref. ft4| . 
For z, z' > 0, we find an identity 



n 1 nX pX pX J /-/ pi J 

CC' / dtt duJ 2v {2tuQJ 2u {2tO =C' 2 / dtt —J 2lJ+1 (2tu()J 2u+1 (2t(') + )-J 2 „(2C) — J 2u+l {2uO , 
Jo Jo Jo Jo u 2 J u 



by a partial integration. The asymptotic formulas (|C4|), together with the identity (|C9|), yield 



S(z,z') ~ 2N 



2 dtt J 2v+ x (2tC) ( 2t( / duJ 2v (2tu() - (2f + 1 



du 



J 2lJ+1 {2tuQ 



[ dtt J 2 „(2iC) (~2tC [ duu J 2l/+1 (2tu(') + (2v + 1) f duJ 2u {2tu£) 
s Jo V Jo Jo 



A2v + i Mf^ l{2<) 

C Jo u 



Partial integrations give rise to the Bessel function equalities 



(2z/+l) / —J 2v+1 (2tu() = -J 2v+1 {2tQ + 2tC / duJ 2u (2tu(), 

Jo u Jo 

(2v + 1) / duJ 2v (2tuC,') = J 2 „(2t(') + 2tC / duuJ 2 u+x(2tuC). 

Jo Jo 



4/ rfUJ 2 , +1 (2tC)J 2 , + i(2K0 + ^4^^ + i(2C)-(2^ + l)^^ ! —hu+x^uQ 
s s Jo u 



Substituting ( Cll ) into ( CIO ) and using the formula (B5), we obtain 
r r 1 

S{z,z) ~ 2N 
By a partial integration, we can rewrite it as 

S(z,z') ~ 27V 



4 / dttJ 2u+1 {2tQJ 2v+1 {2tC) - ' h ^p J C dsJ 2 „ +2 (s) 



(C9) 



(CIO) 



(Cll) 



(C12) 



(C13) 
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This is the asymptotic formula derived by Forrester, Na- 
gao, and Honner E] in a study of parametric random 
matrices and then applied by Akemann and Kanzieper 
|34| to the massive Dirac operator problem. Thus we 
established the equivalence in the case of quadruply de- 
generate masses. We can straightforwardly extend it to 
the formulas with z and/or z' negative. 



[26] 



e-mail addresses: 

nagao@sphinx.phys . sci . osaka-u. ac . jp, 
nishigak@th.phys . titech. ac . jp. 

M.L. Mehta, Random Matrices, 2nd Ed., Academic Press 
(San Diego, 1991). 

T. Guhr, A. Muller-Groeling, and HA. Weidenmiiller, 
Phys. Rep. 299, 189 (1998). 

K. Efetov, Supersymmetry in Disorder and Chaos, Cam- 
bridge Univ. Press (Cambridge, 1997). 
E.V. Shuryak and J.J.M. Verbaarschot, Nucl. Phys. 
A560, 306 (1993). 

J.J.M. Verbaarschot and I. Zahed, Phys. Rev. Lett. 70, 
3852 (1993). 

J.J.M. Verbaarschot, e-print |hep"-ph/990239l . 
P.H. Damgaard, J. Osborn, D. Toublan, and J.J.M. Ver- 
baarschot, Nucl. Phys. B547, 305 (1999). 
D. Toublan and J.J.M. Verbaarschot, Nucl. Phys. B560, 
259 (1999). 

H. Leutwyler and A. Smilga, Phys. Rev. D46, 5607 
(1992). 

A. Smilga and J.J.M. Verbaarschot, Phys. Rev. D51, 829 
(1995). 

J.J.M. Verbaarschot, Phys. Rev. Lett. 72, 2531 (1994). 
J.J.M. Verbaarschot, HA. Weidenmiiller, and M.R. Zirn- 
bauer, Phys. Rep. 129, 367 (1985). 

MA. Halasz and J.J.M. Verbaarschot, Phys. Rev. D52, 
2563 (1995). 

T. Banks and A. Casher, Nucl. Phys. B169, 103 (1980). 
T. Nagao and K. Slevin, J. Math. Phys. 34, 2075 (1993). 
P.J. Forrester, Nucl. Phys. B402, 709 (1993). 
A.V. Andreev, B.D. Simons, and N. Taniguchi, Nucl. 
Phys. B432, 487 (1994). 

K. Slevin and T. Nagao, Phys. Rev. Lett. 70, 635 (1993). 
P.J. Forrester and T.D. Hughes, J. Math. Phys. 35, 6736 

(1994) . 

C.A. Tracy and H. Widom, Commun. Math. Phys. 161, 
289 (1994). 

J.J.M. Verbaarschot, Nucl. Phys. B426, 559 (1994). 
T. Nagao and P.J. Forrester, Nucl. Phys. B435, 401 

(1995) . 

T. Nagao and P.J. Forrester, Nucl. Phys. B509, 561 
(1998). 

T. Guhr and T. Wettig, Nucl. Phys. B506, 589 (1997). 
A.D. Jackson, M.K. §ener, and J.J.M. Verbaarschot, 
Phys. Lett. B387 355, (1996). 

P.H. Damgaard and S.M. Nishigaki, Nucl. Phys. B518, 
495 (1998). 



[27 

[28 

[29 

[30 

[31 

[32 

[33 
[34 
[35 
[36 

[37 

[38 
[39 

[40 
[41 

[42 

[43 

[44 

[45 
[46 



T. Wilke, T. Guhr, and T. Wettig, Phys. Rev. D57, 6486 
(1998). 

J. Jurkiewicz, M. Nowak, and I. Zahed, Nucl. Phys. 
B513, 759 (1998). 

S.M. Nishigaki, P.H. Damgaard, and T. Wettig, Phys. 
Rev. D58, 087704 (1998). 

M.E. Berbenni-Bitsch, S. Meyer, and T. Wettig, Phys. 
Rev. D58, 071502 (1998). 

F. Farchioni, P. de Forcrand, I. H ip, C.B. Lang, and K. 
Splittorff, e-print |hep-lat/9912004[ 



T. Nagao and P.J. Forrester, Nucl. Phys. B530, 742 
(1998); B563, 547 (1999). 

T. Nagao and S.M. Nishigaki, e-print friep-th /000 1137 



G. Akemann and E. Kanzieper, e-print [iep-th/000118g| . 



S. Nishigaki, Phys. Lett. B387, 139 (1996). 

G. Akemann, P.H. Damgaard, U. Magnea, and S. Nishi- 
gaki, Nucl. Phys. B487, 721 (1997). 

E. Kanzieper and V. Freilikher, Philos. Mag. B77, 1161 
(1998). 

H. Widom, J. Stat. Phys. 94, 347 (1999). 

M.K. §ener and J.J.M. Verbaarschot, Phys. Rev. Lett. 
81, 248 (1998). 

F. Dyson, Commun. Math. Phys. 19, 235 (1970) 

M.L. Mehta, Matrix Theory, Les Editions de Physique 
(Paris, 1989). 

G. Mahoux and M.L. Mehta, J. Phys. I (France) 1 (1991) 
1093. 

P.J. Forrester, T. Nagao, and G. Honner, Nucl. Phys. 
B533, 601 (1999). 

T. Nagao and M. Wadati, J. Phys. Soc. Japan 60, 3298 

(1991); 61, 78 (1992); 61, 1910 (1992). 

P.H. Damgaard, Phys. Lett. B424, 322 (1998). 

G. Akemann and P.H. Damgaard, Nucl. Phys. B528, 411 

(1998); Phys. Lett. B432, 390 (1998). 



11 





12 



